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ABSTRACT 

Recent observations by CREAM and ATIC-2 experiments suggest that (1) the spectrum of cosmic 
ray (CR) helium is harder than that of CR proton below the knee 10 15 eV and (2) all CR spectra 
become hard at > 10 11 eV/n. We propose a new picture that higher energy CRs are generated in 
more helium- rich region to explain the hardening (1) without introducing different sources for CR 
helium. The helium to proton ratio at ~ 100 TeV exceeds the Big Bang abundance Y = 0.25 by 
several times, and the different spectrum is not reproduced within the diffusive shock acceleration 
theory. We argue that CRs are produced in the chemically enriched region, such as a superbubble, 
and the outward-decreasing abundance naturally leads to the hard spectrum of CR helium if CRs 
escape from the supernova remnant (SNR) shock in an energy-dependent way. We provide a simple 
analytical spectrum that also fits well the hardening (2) because of the decreasing Mach number in 
the hot superbubble with ~ 10 6 K. Our model predicts hard and concave spectra for heavier CR 
elements. 

Subject headings: acceleration of particles — cosmic rays — shock waves — supernova remnants 



1. INTRODUCTION 

Recently, the Cosmic Ray Energetics And Mass 
(CREAM) has directly observed the CR compositions 
with high statistics in the wide energy range up to 
about 10 14 eV. Interestingly, CREAM shows N P (E) cx 

£-2.66±0.02 for CR proton and jV Hc (£) OC £-2.58±0.02 

for CR helium in the energy region 2.5 x 10 12 eV- 
2.5 x 10 14 eV, that is, the spec trum of CR heli um is 
harder than that of CR proton (|Ahn et all 12010ft . Al- 
though the difference of the spectral index As s» 0.08 
appears small, the implications are of great importance 
as shown below. In addition, the spectral index becomes 
hard by ~ 0.12 for CR proton and by ~ 0.16 for CR he- 
lium at > 2 x 10 11 eV/n because the Alpha Magnet Spec- 
tromete r (AMS) shows NjE ) oc £-2-78±o.oo9 for the CR 
proton (lAlcaraz et al.ll2000aft a nd Nu e (E) cx £-2-74±o.oi 
for CR helium (|Alcaraz et al.ll20 00b) in the low-energy 
range 10 10 eV-10 11 eV. These results have been already 
obtained b y the Advanced Th in Ionization Calorimeter-2 
(ATIC-2) (IPanovet al.ll2009h . 

For CR electrons, the Fermi gamma-ray space tele- 
scope has recently observed the spectrum of CR elec- 
trons in the wide energy rang e from 7 x 10 9 eV to 
10 12 eV (jAckermann et al.|[2010f ). Fermi shows that the 
observed date can be fitted by a power law with spec- 
tral index in the interval 3.03 — 3.13 and the spectral 
hardening at about 10 11 eV, which may have the same 
origi n as that of the CR nuclei. (For other models, see, 
e.g., Kashivama et alJ[2010b iKawana ka et alJl2010t lloka I 



120101 and references therein). Note that we do not dis 
cuss CR positrons in this letter. 

Supernova remnants (SNRs) are thought as the ori- 
gin of the Galactic CRs. The most popular acceleration 
mechani sm at SNRs i s the diffusive shock a cceleration 
(PSA) (lAxford et all [l977 [ IKrvmskvi fl977t iB^lll fl978t 
iBlandford fc Ostriken 11978ft . In fact, Fermi and AG- 



ILE show that middle-age SNRs interacting with molec- 



ular clouds emit gamma-rays (e.g. lAbdo et all [2009; 
iTavani et al.ll2010T ) and the gamma-ray observations sup- 
port that SNRs produce the bulk of Galactic CRs (e.g., 
IQhira et alJl2TTTlt ILi fc Chedl2lHflft . 

According to DSA theory, the spectrum of accelerated 
particles at a shock does not depend on CR elements, 
but depends only on the velocity profile of the shock. 
Thus, naively, recent CR observations seem to show that 
the acceleratio n site of CR helium is different from that 
of CR proton (B lcrmann et alJl2010h . However, in this 
different site scenario, it should be by chance that the ob- 
served ratio of CR helium and proton, N-ft e /N p , at 10 9 eV 
is similar to the cosmic abundance (Y = 0.25). Fur- 
thermore, the difference of the spectral index As ~ 0.08 
means that Nn e /N p at 10 14 eV is about 3 times higher 
than that at 10 9 eV. This enhancement is amazing since 
the mean helium abundance in the universe is virtually 
maintained constant. The stellar nucleosynthesis never 
enhances the mean helium abundance by a factor, which 
is the essential reason that the big bang nucleosynthesis 
is indispensable for the cosmic helium abundance. To 
make the enhancement, we should consider inhomoge- 
neous abundance regions. We show that this leads to 
the different spectrum of CR proton and helium when 
escaping from SNRs. 

In this letter, considering the inhomogeneous abun- 
dance region, we provide a new explanation about (1) 
the different spectrum of CR proton and helium, even 
if CR proton and helium are accelerated simultaneously. 
Our idea uses the fact that CRs escaping from SNRs gen- 
erally have a different spectrum than that of the accel- 
eration site dPtsuskin fc Zirakashv ili 2005; IQhira et all 
I2010L ICaprioli et al.l 12010ft . The runaway CR spectrum 
depends on not only the acceleration spectrum at shocks 
but also the evolution of th e maximum energ y and the 
number of accelerated CRs ([Ohira et alj [2010). We also 
suggest that (2) the spectral hardening of CRs is caused 
by the decreasing Mach number in the high tempera- 



2 



Ohira and loka 



Runaway CR 



(s+a//3) 




Energy 



Fig. 1. — Schematic picture of the runaway CR spectrum. The 
solid, dashed and dotted lines show the runaway CR spectrum, 
CR spectrum inside an SNR at an early epoch and CR spectrum 
inside the SNR at a later epoch, respectively. The solid line of the 
runaway CR spectrum represents the equation J3J. A variable \ 
(e.g., the shock radius) describes the SNR evolution. 

ture medium. Both the inhomogeneous abundance and 
the high temperature can be realized in the superbubblcs 
with multiple supernovae. Our conclusions are summa- 
rized as follows. 

• Runway CR spectra depend on not only CR spec- 
tra inside the SNR but also the evolution of the 
maximum energy and the number of accelerated 
CRs. Therefore, taking account of the inhomoge- 
neous abundance region, runaway CR spectra of 
different CR elements have different spectra (sec- 
tion 2 and 3.1). 

• Our model is in excellent agreement with observed 
spectra of CR proton and helium. Harder spec- 
trum of CR helium is due to the enhancement of 
the helium abundance around the explosion center. 
On the other hand, the concave spectra of all CR 
elements are due to the decreasing Mach number in 
the hot gas with ~ 10 6 K. The concave spectra may 
be also produced by the CR nonlinear effect, the en- 
ergy dependent effects on the accelerated CRs (on 
a or ff), the propagation effect (7), and/or multi 
components with different spectral indices (section 
3.2 and 4). 

• Within the single component scenario, the hard 
helium spectrum suggests that the origin of the 
Galactic CR is SNRs in superbubbles, although 
we are not excluding the multi component scenario 
(section 5). 

• Our model predicts that heavier (at least volatile) 
CR elements also have harder spectra than that of 
CR proton and have concave spectra (section 5). 

2. RUNAWAY CR SPECTRUM 

In this section, we br iefly review the run away CR spec- 
trum (see Appendix of lOhira et al. (2010)). We here use 
a variable x (for example the shock radius or the SNR 
age) to describe the evolution of an SNR. Let .Fsnr(XjP) 
and p m ax(x) be the CR momentum spectrum [(eV/c) -1 ] 
and the maximum four momentum of CR inside the SNR 
at a certain epoch labeled by \: respectively. CRs escape 
in order, from the maximum energy CR because the dif- 
fusion length of high-energy CRs is larger than that of 



low-energy CRs. Then, the number of runaway CRs be- 
tween x and X + dx is 



Fsnr(x,P* 



dp n 



-dx 



(1) 



which corresponds to the number of runaway CRs be- 
tween p = p m ax(x) and P = Pmax(x) + dp, F csc (p)dp. 
Hence, F csc (p) is 



F esc (p) = F SNR (p m i x (p),p) 



(2) 



where p max (p) is the inverse function of p max (x)- Assum- 
ing F SNR (x,p) oc x P P~ s and p max (x) oc x~", we obtain 
the runaway CR spectrum as 

Fc S c(p)°cp-( s +£) , (3) 

where a and ft are parameters to describe the evolution 
of maximum energy and the number of accelerated CRs, 
respectively. (We use a ~ 6.5 and (3 ~ 1.5 later.) There- 
fore, the runaway CR spectrum F csc is different from that 
in the SNR Fsnr oc p~ s . Figure Q] shows the schematic 
picture of the runaway CR spectrum. In this Letter, we 
use the shock radius, i? s h, as x- 

The evolution of the maximum energy of CRs at 
the SNR has not been understood. This strongly de- 
pends on the evolution of the magnetic fi eld around 
the shock (e.g. iPtsuskin fc Zirakashvilill2003l ). Although 
some magnetic field amplifications have been proposed 



2007; 



Lucek fc B ell 2000; Bell 20Ql iGiacalone fc Jokipiil 



Ohira ct al. 2009b) and investigated by simulations 



Niemiec et al 1120081 iRiauelme fc Spitkovskvi [2009; 



lOhira et al.ll2009aHGargate et al.ll2010D . the evolution of 
the magnetic field has not been completely understood 
yet. Here we assume that CRs with the knee energy es- 
cape at R = i?sedov, where i?sodov is the shock radius 
at the beginning of the Sedov phase. Furthermore, we 
use the phe nomenological approach with the pow er-law 
dependence (jGabici et all 120091: lOhira et alj|2010f >. 



(i?sh) 



V-Rscdo- 



(4) 



where pknoc = 10 15 ' 5 eV/c is the four momentum of the 
knee energy. Note that a does not depend on the CR 
composition because the evolution of the maximum en- 
ergy depends only on the evolution of the magnetic field 
and the shock velocity. 

The evolution of the number of CRs inside the SNR 
has not been also un derstood. This d epends on the 
injection mechanism (jOhira et al.l [20 H ) :■ and the den- 
sity profile around the SNR. We here ad opt the ther- 
mal leakage model (jMalkov fc Volld I1995D as an injec- 
tion model. For the total density profile, ptot{Rsh) ~ 
m p (n p (i? s h) + 4nHc(-Rsh)) where n p and nn c are the 
number density of proton and helium and m p is the pro- 
ton mass, the shock velocity of the Sedov phase is 



Wsh(-Rsh) OC Ptot(-Rsh) 2 R : 



sh 



(5) 

In the thermal leakage model, the injection momentum 
of element i is proportional to the shock velocity, p- m j,i oc 
u s h, and the number density of CR with momentum Pi n j,i 
is proportional to the density, Pi niii fi{Pinj,i) oc rii(R sh ), 
where fi is the distribution function of CR element i. 
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Fig. 2. — Schematic picture of the formation of the different spec- 
trum. The solid and dashed line show the proton density and the 
helium density, respectively. The dotted lines show the shock front. 
In early phase, high-energy CR proton and CR helium escape, and 
in late phase, low-energy CR proton and CR helium escape. The 
ratio of CR helium to CR proton increases with the CR energy. 

Hence, the number of CR element i with a reference mo- 
mentum p = m p c, i*sNR,i(-Rsh,«T-pc) is 



FsNR,i(Rsh, m p c) oc R^ h fi(m p c) 



ccni(R sh )p tot (R sh ) 2 ° w i? s " 



3(3- 



where fi(p)p 2 oc p S1 



and si ow is the spectral index 
in the non-relativistic energy region. For the nonlinear 
DSA, the spectral index in the non-relativistic energy 
regio n is different from that in the relativistic energy re- 
gion (Bcrczhko & Ellisonl fl999[) . To understand the es- 
sential feature of the runaway CR spectrum, we here con- 
sider only the test-particle DSA, that is, s\ ow — s. Be- 
cause ni(i? s h)/Ctot(^?sh) _5_ is not always a single power- 
law form, the evolution of the number of accelerated CRs 
can not be always described by a constant (3. 

3. BASIC IDEA 

3.1. Different spectrum of CR proton and helium 

According to the test particle DSA theory, the index s 
of relativistic CR energy spectrum depends only on the 
velocity jump at the shock, 



ui + 2u 2 
u\ - u 2 



M 2 



1 



M 2 -l 



(7) 



where we use the Rankinc-Hugoniot relation at the sec- 
ond equation and M is the Mach number. Then, the 
index of the runaway CR spectrum, s osc , is 



ft 

s H 

a 



(8) 



in equation ([3]). Therefore, if j3/ a (in particular f3, the in- 
dex for the accelerated CR number evolution) is different, 
the runaway CR spectrum is different between the CR 
compositions. This is our main idea to explain the he- 
lium hardening observed by CREAM and ATIC-2. From 
equation (|6]), (3 depends on the ambient number density 
m. Therefore, different density profiles make different 



runaway CR spectra (See Section @] for more details). 
Figure [2] shows the schematic picture of our idea. 

3.2. Spectral hardening of all CRs at the same energy 
per nucleon 

In this subsection, we discuss the spectral hardening of 
the observed CRs. The Galactic CR spectrum observed 
at the Earth, -Fobs, is obtained by the simple leaky box 
model 



Fobs oc F csc (p)/D(p) oc F esc (p)p 7 



(9) 



where D(p) oc p 7 is the diffusion coefficient (e.g. 
IStrong et aLl[200l . 
spectrum is 



Hence, the index of the observed 



^obs 



a 



(10) 



The deviation from a single power law means that at 
least one of s, a, j3, and 7 has an energy dependence 
or that the origin of low energy CRs below 10 11 eV is 
different from that of high energy CRs above 10 11 eV. 
Although the multi component scenario may be the case 
because there are many types of SNRs, we discuss the 
single component scenario in this letter. 

Firstly, we discuss the energy dependence of s. From 
equation (0), s depends on the shock radius because the 
Mach number M decreases with the shock radius. Then 
we can expect the spectral harding of all CR compo- 
sitions at the same rigidity cp/Ze, that is, at approxi- 
mately the same energy per nucleon. From equation ([5]) , 
the Mach number is 



M w 10 d 



Ptot(-Rsh) 
Ptot(-Rsedov) 



T 



10 4 K 



/ Rsh 



\Rs- 



cdov 



(11) 

where T is the surrounding temperature and we assume 
that the ejecta mass and the energy of supernova explo- 
sion are 1 M and 10 51 erg, respectively. From equa- 
tions (HJ), ([7]) and pT|) . we can obtain s as a function of 
p (see §11]). 

Alternatively the spectral hardening ca n be also inter- 
preted as the CR nonlinear effect (e.g., IDrurv fc Volkl 
HMD: IMalkov fc Drurvl l200l . This issue will be ad- 
dressed in the future work. 

Next, we discuss the energy dependence of (3 that is 
the parameter to describe the evolution of the number 
of accelerated CRs. In Section [3~T1 we consider different 
power-law forms for n p (i? s h) and nn e (R s ] 1 ) to make the 
different spectrum of the CR proton and helium. There- 
fore, jOtot(-Rsh) ~ rn p [n p (R sh ) + Ann c (R sh )] is not a single 
power law form, and (3 has an energy dependence (see 

m- 

The energy depende nce of 7 will be s oon precisely de- 
termined by AMS-02 (|Pato et al.ll2010[ ). We do not dis- 
cuss the energy dependence of a because the complete 
physics of the CR escape and magnetic turbulence is be- 
yond the scope of this Letter. 

4. COMPARISON OF OUR MODEL WITH OBSERVATIONS 

In this section, specifying model parameters, we calcu- 
late the Galactic CR spectrum. For simplicity, we here 
assume the number densities of proton and helium as 
follows, 

np(Rsh) = np,o 
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Fig. 3. — Comparis on of our model (so lid and dashed line) 
with AMS (triang le) (Alcaraz ct al. 2000 aB), ATIC-2 ( square) 
UPanov et al.ll2009D and CREAM (circle) < Aim et alj|2uTO) obser- 
vations for Galactic CRs. Filled symbols and the solid line show 
CR proton. Open symbols and the dashed line show CR helium. 
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Fig. 4. — The evolution of the irieixirQurQ energy .c/max 

jZ (solid 

line) and the spectral i ndex s of test particle DSA (dashed line) in 
equations Q, J7) and JTTJ as functions of R a ^. 



n He (i?sh) = C«p, 



/ Rsh 



V -Rscdo 



(12) 



where n P) o is the number density of proton at i? s h = 
-Rsedov, and C n p,o is the normalization factor of the he- 
lium density. We set ( = l0 6 - 5( . s / a )^ 1 so that the helium 
abundance is that of the solar abundance, nn c /n p = 0.1 
(i.e., Y w 0.25), when cp max = Z GeV with equation 
((U). Note that the power-law dependence is a first step 
approximation for the mean value. Then, from equa- 
tions ((2|), (j4j, ©, ©, observed spectra of CR proton 
and helium are 



-^obs,p — -^p,knee 



1 + C (p/p 



knee J 



i + C 



Pkn 



3{3-s(p)} 



+ "'] 



, (13) 



-^obs, He — ^-fp,knee 



1 + C(p/ Z Pknee)' 



i + C 



Zp 



kne 



-U(p)+ 3<3 - s £ )} - 25 +-r] 



(14) 



where -Fp^kncc and e-Fp,knee are normalization factors of 
CR proton and helium and Z — 2 for helium, and s(p) is 
obtained from equations ((4]), ([7]) and ([11]) . In this model, 
all parameters are a, 7, S, e, T, -F^knee- 

Figure [3] shows the comparison of our model with 
observations. We take into account the solar modula- 
tion effects with the mod ulation potential $ = 450 MV 
(jGleeson fc Axfordl U968h . Our model is in excellent 
agreement with the observed spectra, with a = 6.5, 
7 = 0.43, S = 0.715, e = 0.31, T = 10 6 K. The dif- 
ferent spectra of CR proton and helium originate from 
the different density profiles in equations (fT"2"l). Figured] 
shows the evolution of the maximum energy of CRs and 
the spectral index of CRs inside the SNR. In the early 
phase, the spectral index s is 2 and after then, the spec- 
tral index decreases with the shock radius because the 
Mach number decreases with shock radius. The change 
of spectral index s is about 0.1 which is almost the same 
as the observed hardening. The observed hardening is 
not the result of the change of the injection history, /?, 
but the result of spectral change of CRs inside the SNR. 
The high temperature T ~ 10 6 K is necessary for the 
spectral hardening As ~ 0.1. 

In addition, our model also makes a concave spectrum 
of CR electrons as observed (|Ackermann et all 120101) . 
However, the evolution of injection efficiency of CR elec- 
trons has not been understood well. So we need further 
studies to discuss the CR electron spectrum in detail. 

5. DISCUSSION 

To make the different spectrum, our model requires 
that the helium abundance around the explosion cen- 
ter is higher than that of the solar abun dance. SNRs 
in sup erbubbles are one of candidates. iHigdon et al.l 
(1998) show that supernova ejecta can dominate the 
superbubble mass within a core radius of one third of 
the superbubble radius. In the stellar wind and the 
supernova explosion, the stellar hydrogen envelope has 
lower density and higher velocity than that of helium. 
Then we expect that the helium fraction in the cen- 
ter of superbubbles is higher than that in the outer re- 
gion. Furthermore, to make the concave spectrum, our 
model requires an ambient medium with high temper- 
ature, T = 10 6 K. This is also consistent with su- 
perbubbles. According to the CR composition study, 
SNRs in superbubbles hav e been considered as the ori- 
gin of Galact i c CR s (e.g.. [Lingcnfcl ter fc Higdonl 120071 
lOgliore et al.l 120091 ) . Particle accelerations in super- 
bubbles have been also investigated by intensive stud- 
ies ("e.g.. iBvkov fc Fleishman! Il992t iParizot etall 12004 
iDar fc De Ruiula I 120081: iFerrand fc Marcowithl 120101) . 
We here considered a spherically symmetric system. The 
off-center effects may be important for the initial phase 
and thereby for the high energy spectrum, because the 
shock radius at the beginning of the Sedov phase -Rsodov 
is about 20 pc which is c omparable to the ty pical size 
of OB association, 35 pc (|Parizot et al.|[20CM) . and the 
shock radius i? s h is about 200 pc at the end of the Sedov 
phase. This is an interesting future problem. 

Note that the spectral hardening can be also made by 
the nonlinear model, the energy dependence of the CR 
diffusion coefficient and/or multi components with dif- 
ferent spectral indices. So the high temperature may not 
be absolutely necessary. The stellar wind of red giants is 
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one of candidates for the cold and helium rich ambient. 
Still, the dominant core -collapse supernovae is type II 
(e.g.. lSmartt et aLll2009D which has no helium rich wind, 
so that the superbubble scenario looks more likely as the 
origin of the Galactic CRs above 10 11 eV. For the CRs 
below 10 11 eV, the spectral difference between CR pro- 
ton and helium may be cause d by the solar mo dulation 
and the inelastic interactions (|Putze et al.ll2010h . 

The spatial variation of the helium ionization degree 
can also change the injection history. The injection effi- 
ciency of the large rigidity is thought to be higher than 
that of low rigidity since particles with large rigidity can 
easily penetrate through the shock front from the down- 
stream region. If the ionization degree increases with the 
SNR radius, the CR helium spectrum becomes harder 
than the CR proton one, /3ne < ftp- However, the rigidity 
dependence of the injection efficiency has not been under- 
stood completely. Moreover, the inj ection from neutral 
parti cles should also be understood (jOhira "et~aH 2009b. 

According to our model, CR spectra of heavier volatile 
elements than helium is also harder than that of proton. 
Low-energy CRs of refractory elements are thought to re- 
sult from suprathermal injectio n by sputtering off p reac- 
celerated, high-velocity grains (jEllision et al.1 11997ft . To 
be accelerated to the relativistic energy, the refractory 
elements should be sputtered because the grains can not 
be accelerated to the relativistic energy. The SNR shock 
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velocity is not fast enough to accelerate refractory ele- 
ments to the knee energy when the refractory elements 
are injected because the sputtering time scale is too long. 
Therefore, refractory CRs around the knee energy should 
be injected by the standard manner similar to volatile 
CRs. In this case, the refractory CRs also have harder 
spectra than protons, although we need further studies 
of the injection of refractory CRs at the knee energy. 

If CRs trapped inside the SNR and released at the end 
of t he SNR's life outnum ber runaway CRs (see figure 
3 in iCaprioli et al.l (|2010f ) ) , our scenario does not work 
for producing hard and concave spectra. In our model 
with a ~ 6.5 in Eq. trapped CRs have energy below 
1 GeV when they are released, that is, p max J$ Zm p c 
when i? s h > 10i?sodov, and are not relevant for our in- 
terest. Higher energy CRs escape from the SNR even 
after advected to the downstream since the CR diffusion 
is faster than the expansion of the SNR. Our case is sim- 
ilar to the right figure 7 in ICaprioli et al.l (|2010() where 
trapped CRs are released below 100 GeV. The energy 
boundary between trapped CRs and runaway CRs de- 
pends on the evolution of the maximum energy (a) . 
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